Optical conductivity of a metal-insulator transition 
for the Anderson-Hubbard model in 3 dimensions away from 1/2 filling 



X. Chen and R. J. Gooding 
Department of Physics, Queen's University, Kingston, ON K7L 3N6, Canada 

(Dated: August 14, 2009) 

The Anderson-Hubbard model is considered to be the least complicated model using lattice 
fermions with which one can hope to study the physics of transition metal oxides with spatial 
disorder. We have completed a numerical investigation of this model for three-dimensional simple 
cubic lattices using a real-space self-consistent Hartree-Fock decoupling approximation for the Hub- 
bard interaction. In this formulation we treat the spatial disorder exactly, and therefore we account 
for effects arising from localization physics. We have examined the model for electronic densities 
well away 1/2 filling, thereby avoiding the physics of a Mott insulator. Several recent studies have 
made clear that the combined effects of electronic interactions and spatial disorder can give rise to 
a suppression of the electronic density of states, and a subsequent metal-insulator transition can 
occur. We supplement such studies by calculating the ac conductivity for such systems. Our nu- 
merical results show that weak interactions enhance the density of states at the Fermi level and 
the low-frequency conductivity, there are no local magnetic moments, and the ac conductivity is 
Drude-like. However, with a large enough disorder strength and larger interactions the density of 
states at the Fermi level and the low-frequency conductivity are both suppressed, the conductivity 
becomes non- Drude-like, and these phenomena are accompanied by the presence of local magnetic 
moments. The low-frequency conductivity changes from a a — Ode ~ uj 1 ^ 2 behaviour in the metallic 
phase, to a a ~ u> 2 behaviour in the nonmetallic regime. For intermediate disorder at 1/4 electronic 
filling, a metal-to-insulator transition is predicted to take place at a critical U /B « 0.75 (U being 
the Hubbard interaction strength and B the electronic band width). Our numerical results show 
that the formation of magnetic moments is essential to the suppression of the density of states at 
the Fermi level, and therefore essential to the metal-insulator transition. At weaker disorder a small 
lessening of the density of states at the Fermi level occurs, but screening suppresses the spatial 
disorder and with increasing interactions no metal-insulator transition is found. 



I. INTRODUCTION 

Electron-electron interactions are important in un- 
derstanding the properties of many transition metal 
oxides. 1 ^ Novel ordered phases are found in this class 
of compounds: e.g., high-T c superconductors, d-electron 
heavy fermions, and quantum magnets. The inclusion of 
disorder into such materials, and into the models of such 
systems, adds another level of complexity. In this report 
we focus on metal-to-insulator transitions that are con- 
trolled by the combined effects of both interactions and 
disorder. In part, this work is part of an effort to under- 
stand several experimental results, including the proper- 
ties of the weakly doped cuprates,-' 4 for which disorder 
effects arc known to be important. 

Initially, we were motivated to conduct such theoreti- 
cal work to better understand this transition in the ma- 
terial LiALjTia— The ordered and undoped mate- 
rial LiTi204 undergoes a transition to a superconduct- 
ing phased around T c w 12 K, and has been proposed 
to be related to the family of high-T c superconductors^ 
If such a conjecture is correct, there must be reasonably 
large electron-electron interactions in this system. The 
Al-doped system undergoes a metal-insulator transition 
for y « 0.33, and recent theoretical work by one of us and 
co-workers have suggested that (i) the effects of disorder 
alone cannot lead to such a transition^; (ii) inclusion of 
electron-electron interactions in the form of an on-site 
Hubbard energy can lead to such a transition*^ This oc- 



curs by the suppression of the density of states at the 
Fermi energy to near zero. The results presented in this 
report are a natural continuation of such work. 

In addition, other systems display such transition. 
Sarma et aliii have investigated the electronic structure 
of LaNii_ :E Mn :E 03. It is known that for a critical con- 
centration of i c ~ 0.1 a metal-to-insulator transition is 
found. The point-contact tunnelling conductance spec- 
tra, which provides information on the density of states, 
show a suppression at the Fermi energy in the form of 
a downward cusp for both x = 0.05 and 0.1. The con- 
ductance is found to be proportional to VV, where V 
is the bias voltage. This behaviour was predicted by 
Altshuler and AronoA*^, and is similar to the results 
that we encountered in studies of LiAL y Ti 2 -_j,04i 1 £ Kim 
et al»i2 have studied the transport and optical proper- 
ties of SrTii„ x Ru x 03. When x decreases from 1 to 0, 
SrTii-^Ru^Os evolves from a correlated metal, SrRuOa 
(x = 1), to a band insulator, SrTi03 (x — 0). Depend- 
ing on x, there are six types of electronic states. The 
metal-to-insulator transition takes place at x c ~ 0.5. As 
x decreases from one, the concentration of conduction 
electron decreases, and the effective interaction strength 
increases. 

These systems are very complicated, having several 
important d orbitals for different conducting sites, and 
many different interaction, disorder and hopping energy 
parameters. The theoretical considered below are a min- 
imalist's approach to such interesting behaviour. 
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A. Anderson-Hubbard Model: 

A simplified model that describes interacting electrons 
moving on a spatially disordered lattice is the so-called 
Anderson-Hubbard model. Its Hamiltonian is given by 

(ij),a i,v i 

(1) 

Electron annihilation/creation operators for site j and 
spin a =T, J. are represented by c JCT and ct , respectively. 
The spatially inhomogeneous environment in which the 
electrons move is accounted for by on-site energies, £j, 
which are random. Often, these are selected to be chosen 
from a uniform distribution, e% £ [— and therefore 

is the energy scale characterizing the strength of the 
disorder. Interactions between electrons are accounted 
for by the intra-site Hubbard interaction term, charac- 
terized by the energy scale U. The near- neighbour hop- 
ping frequency is denoted by t, and all other symbols have 
their usual meaning. It is hoped that this model can cap- 
ture some of the essential physics of the metal-insulator 
transitions 2 of disordered transition-metal oxides. 

The electronic properties of systems described by this 
model Hamiltonian are indeed complicated, as can be 
understood from the following reasons. When disorder is 
sufficiently strong it can lead to the localization of elec- 
trons, and favours large electronic occupancies on sites 
with low on-site energies. Therefore, disorder generates 
both localization effects and an inhomogeneous distribu- 
tion of electronic charge. There exists a critical disorder, 
(W/t) c , beyond which all states are localized and the sys- 
tem becomes an Anderson insulator.— In a spatially uni- 
form half-filled system, electron-electron interactions also 
lead to the effective localization of electrons; however, in 
contrast to disorder, interactions favour single occupancy 
of electrons on all sites. In general, at a critical interac- 
tion, (U/t) c , the system undergoes a transition from a 
metal to a Mott insulator. However, the Mott transition 
may not take place when the electronic concentration is 
away from half-filling. 

Due to the randomness of the on-site energies, it 
is usual to address this model using numerical tech- 
niques. Also, except for some small clusters that 
can be solved exactly, large lattices have to be solved 
with the help of approximate schemes. This diffi- 
culty notwithstanding, a small number of "exact" nu- 
merical results are available, coming from both ex- 
act diagonalization (ED^ii&iLi&i^ and Monte Carlo 
(MC)^.^^^^!^. studies. Because of their 
relevance to our work, we mention one aspect of the re- 
sults in two of these papers. That is, some ED calcula- 
tions showed the presence of a suppression of the DOS at 
the Fermi energy, in both one 2 ^ and two 2 ^ dimensions. 
These results are consistent with the HF papers men- 
tioned below. Therefore, this provides partial verification 
of the results based on the HF method, the method used 
in the remainder of this paper. 



B. Discussion of Previous HF Results: 

The Hartree-Fock (hereafter HF) method is well 
known, and its application to disordered systems is exten- 
sive. The accuracy of the results was recently critiqued 
by us and coworkers where it was shown that provided 
one allowed for sufficient magnetic degrees of freedom the 
energies and charge densities of (small) exactly solvable 
systems agreed well with those obtained from HF. How- 
ever, since the spin correlations are essentially those of 
pairs of classical spins of variable directions and lengths, 
and therefore do not include quantum fluctuations, HF 
is less successful at capturing the correct spin-spin corre- 
lations. 

As mentioned, the HF approximation has been ap- 
plied in many studies of the Anderson-Hubbard and 
related models. Some examples are (i) a study of 
a two-dimensional point-defect model 3 - which repre- 
sents the acceptors and donors in the high-T c cuprate 
La2_ 2; Sr a ,Cu04 through determining the magnetic phase 
diagram, (ii) a proposal of a novel inhomogeneous metal- 
lic phase in two dimension a 32 i 33 ' 34 i 35 as a combined ef- 
fect of disorder and electronic interactions when their 
strengths are comparable, and (iii) a detailed study of 
the three-dimensional Anderson-Hubbard model, deter- 
mining both magnetic and electric phase diagrams at 
half-filling, 3 ^ 

In Refs i 36 ' 37 , Tusch and Logan have focused on 1/2- 
filling and some other fillings lower than 1/8. At 1/2- 
filling with a disorder strength of W/t — 5, the density of 
states (hereafter denoted by DOS) shows a suppression 
for both U/t = 6 and 9. Using the inverse participa- 
tion ratio (IPR) technique, the system is determined to 
be metallic at U/t = 6 and insulating at U/t = 9. The 
IPR is compared with a threshold mean IPR (obtained in 
the noninteracting limit) that scales with system size to 
determine whether the system is metallic or insulating. 
However, the effect of the interactions on the threshold 
mean IPR is unknown, as are the effects the unusual sta- 
tistical distribution of the IPR, the latter having been 
discussed by Mirlin. 38 In part to circumvent such prob- 
lems, here we focus on the characterization of the metal- 
licity of a system via its optical conductivity. 

In the HF treatment of LiAl y Ti2_ y 04 by one of us 
and coworkers^ a metal-to-insulator transition, again 
found from an examination of the behaviour of the IPR, 
was found. In that work the local magnetic moments 
were restricted to lie along the z-axis, geometrically frus- 
trated corner-sharing tetrahedral lattices were studied, 
and a quantum site percolation model of disorder was 
used. Open questions concern whether the suppression 
of the DOS at the Fermi level depends on the lattice type 
and the disorder model. To answer these questions, we 
choose to study three-dimensional simple cubic lattices, 
which arc unfrustrated, and consider a uniform box dis- 
tribution type of disorder. We also allow local magnetic 
moments to develop in the xz-plane, which increases the 
spin degree of freedom, believed to be important in some 
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circumstances^ 

We also mention the effective-field theory analysis 
of local moment formation in disordered interacting 
Fermi liquids by Milovanovic, Sachdev, and Bhatt. 39 
The potential importance of such moments in the metal- 
insulator transition was left as an outstanding question, 
and will be one of the main aspects addressed in this 
paper. 

C. Summary of New Results: 

In this paper, we apply the real-space self-consistent 
HF method to electrons moving on simple cubic lattices 
at 1/4-filling, with various strengths of interaction and 
disorder. Here we report our results from calculations of 
the DOS and ac conductivity. 

When examining systems with a disorder strength of 
W/t = 6, we find that the DOS and the low-frequency 
conductivity are enhanced by a weak interaction (U/t < 
3), and for this range of interactions there are no mag- 
netic moments in the system, and the ac conductivity 
is Drude-like. With a stronger interaction (U/t > 4), a 
suppression of the DOS at the Fermi level and qualita- 
tive changes (non-Drude-like) in the low-frequency con- 
ductivity are found. We find that concomitant with these 
changes in behaviour is the appearance of local magnetic 
moments in the system, although no evidence of mag- 
netic ordering is present. For this disorder strength and 
electron concentration a metal-to-insulator transition is 
likely to take place at a critical U c /t w 9; that is, roughly 
3/4 of the noninteracting bandwidth. We have also ex- 
amined the weaker disorder strength of W/t — 2, and 
although one finds a small suppression of the DOS with 
increasing Hubbard interaction, for no value of U/t do 
we find a metal-to-insulator transition. 



II. REAL-SPACE SELF-CONSISTENT 
HARTREE-FOCK APPROXIMATION 

For ordered systems, one may transform the Hamilto- 
nian to wave vector space, and then expand the many- 
particle wave function in a complete set of Bloch wave 
functions in the corresponding periodic potential. How- 
ever, in disordered systems there is no translational sym- 
metry, and working in wave vector space does not sim- 
plify the problem. Therefore, we will work in real space. 
In the real-space formulation of HF theory (e.g., see page 
349 of Ref4£) the disorder is treated exactly, and there- 
fore this approximation allows for us to include in our 
calculations the effects of localization. In the real-space 
formulation of HF theory the local Hubbard interaction 
term is replaced by 

w (n a )n 2T + (n iT )n;x - (2) 



Here, the terms that are proportional to fluctuations 
about the mean values squared are ignored. Also, we 
do not consider the expectation values that arise from 
superconducting correlations. 

Substituting Eq. J2]) into the Anderson-Hubbard 
Hamiltonian in Eq. (JTJ), one finds the HF effective Hamil- 
tonian given by 

w eff =-*E fe ->, CT (4s CT +4o 

+ E,>e?aU„ - uY,i(K^ii + h M&r)> ( 3 ) 

where 

e^ = ei + Un lS n ia = (n irT ) (4) 

K = (c^n) = K = <c|^T> = <^ r ) (5) 

Here, m a — (hi a ) is the expectation value of the number 

rr 

of electrons with spin a on site i; ef^ is the spin depen- 
dent effective on-site energy at site i, and it is the sum of 
the original on-site energy plus the Hubbard U times the 
number of electrons with the opposite spin on that site. 
Therefore, in the absence of the last term one includes 
the effects of the Hubbard interaction via an effective 
shift of the on-site energies. The last term in Eq. © 
does include the hf effective local fields, and these guar- 
antee that the effective Hamiltonian is invariant under 
rotations in spin space. In the self-consistent formula- 
tion of HF theory, one must ensure that the solutions of 
the effective one-electron Hamiltonian given in Eq. ([3]) 
lead to local spin-resolved charge densities and effective 
local fields that satisfy Eqs. (|4I5|) when the expectation 
values are taken with respect to the HF ground state 
wave function. 

A. Numerical Approach: 

As mentioned above, one is required to solve the ef- 
fective Hamiltonian self-consistently. In a system with 
N s sites, one has AN S variational parameters that have 
to be determined numerically so that the ground-state 
energy is minimized. (In addition, one includes a chem- 
ical potential to fix the electronic density at whatever 
concentration is desired. In our results below we focus 
on 1/4 filling.) To solve the effective Hamiltonian self- 
consistently, one begins with a random initial guess for 
the AN S parameters that satisfies the constraint of fixing 
the total number of electrons, and then iterates to self 
consistency. One needs to repeat the above procedure 
with many other initial guesses of the variational param- 
eters and below we examine the self-consistent states hav- 
ing the lowest energy HF energies, -Ejjjr- 

After self-consistency is achieved for the HF solutions, 
we obtain effective single-electron energies, {E n }, and 
the DOS is then given by 

DOS(E) = ^S(E-E n ). (6) 

n 
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In our results shown below, we broaden each delta func- where Xq is where fix) takes maximum value, and d is 
tion using a Gaussian function given by the standard deviation (or the broadening width). 

/(*) = 1 e-t'--' (7) 
ay 2n 

I 



To calculate the ac conductivity, one first obtains the imaginary part of the current-current correlation function 
using the Kubo formula: 

ImA xx (q = 0,u;) = -^ J2 y n ' m (f(tm)-f(e n ))6(Lo-(e m -e n )), (8) 



where 



i j,<T,er' 



Then, the real part of the ac conductivity i 



,41 



]mA xx (q=0,u) 

C7l(w) = . (10) 



To calculate the ac conductivity numerically, we broaden each delta function in Eq. ([5]) using a Gaussian function 
given by Eq. ([7]). In our results below we state the broadening for each quantity. 

I 



III. RESULTS 
A. Disorder Strength: 

While the Anderson-Hubbard model is often used to 
study disordered electronic systems with strong electron- 
electron interactions, it is not always clear how the pa- 
rameter space studied for such models relates to real 
materials, particularly with regards to the strength of 
the disorder potential. Here we focus on the appropriate 
range of disorder strengths. 

When characterizing the strength of the potential en- 
ergy it is usual to compare U to the noninteracting band- 
width for electrons moving on an ordered lattice. For a 
3d simple cubic lattice that is B = 12t. Weak interac- 
tions correspond to U/B <C 1, intermediate to U/B ~ 1, 
and strong to U/B ^> 1. This leads to the question, is 
it permissible to use the same characterization for the 
disorder strength? 

We have examined the U = Anderson model (at 1/4 
filling - the same as the electronic concentration studied 
throughout this paper) using a box distribution disor- 
der - all on-site energies between — W/2 and +W/2 are 
equally probable. We have then found the ac conductiv- 
ity of this system for various W/t, and here we discuss 
the results corresponding to W/t =1,2 and 6. If one fits 
our conductivity data to the real conductivity of a Drude 
model, the latter given by 

ai &) = i rr, ^2 ' ( n ) 

1 + [IjJtY 



TABLE I: Conducting properties of noninteracting electrons 
moving on a 3d simple cubic lattice with a box distribution of 
disorder of strength W, for an electron filling of 1/4. The adc 
and t are found from fits of the ac conductivity to Eq. (|11[) . 



W/t 


a dc {iTe 2 t 2 } 


r [i- 1 ] 


1 


15. 


11. 


2 


2.3 


3.6 


6 


0.26 


0.45 



dc conductivity shown in Table Q] The difference be- 
tween the ijj — > extrapolations quantifying the dc con- 
ductivities are striking, in that for W/t — 1 we find 
adc ~ 157re 2 £ 2 , whereas for W/t — 6 we find adc ~ 
0.267re 2 i 2 , implying that the conductivity of these sys- 
tems differs by almost a factor of 60. Clearly, the non- 
interacting W/t = 6 (or equivalently W/B = 1/2) sys- 
tem corresponds to a "bad metal" . (However, while the 
W/t — 6 system is a very poor conductor, for this disor- 
der strength we are not approaching a 1/4-filled Ander- 
son insulator, and, in fact, our calculations and existing 
datai^ allow us to determine that ep — nc ~ 4t, with [ic 
being the location of the mobility edge.) 

Further justification for this claims follows from a sec- 
ond way of characterizing the "strength" of the disorder, 
and corresponds to evaluating the elastic mean free path, 
l ei in the Born approximation. One finds 42 

l - « - (-Y (12) 
a z 2 \WJ ' 1 ] 



one obtains estimates of the relaxation time and the 



where a is the lattice constant and z is the coordi- 
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nation number for a particular lattice. Therefore, for 
W/t = 1, 2, and 6 one has l e /a ~ 50, 13, and 1.4, re- 
spectively. While one does not expect the Born approx- 
imation to be accurate quantitatively, especially when it 
predicts £ e rs a, from these results it follows that for 
nonintcracting electrons moving on a three dimensional 
simple cubic lattice with W/B = 1/2 one has a very 
short mean free path. So, while it might be conventional 
to refer to a U/B = 1/2 correlated electronic systems 
as being intermediate coupling, for a disorder strength 
of W/B = 1/2 the system remains metallic but the ef- 
fect of the disorder potential is indeed very large. As we 
see in the results below, it is only for such large disorder 
strengths, which by itself leads to substantial localiza- 
tion effects, that we find a disorder+correlation-driven 
metal-insulator transition. 

(Also, note that in our previous study of lithium 
titanate^ we employed a quantum site percolation 
model, which formally corresponds to setting all A site 
energies to be zero while those on the B sites are infinite 
- that is, the B sites are removed from the conduction 
path. With such a large disorder potential (of infinite 
strength for a binary alloy model of disorder) we did ob- 
tain a metal-insulator transition^) 



B. Variation of the Density of States: 

As mentioned in the introduction, one of the motiva- 
tions for this study was to better understand the con- 
ditions necessary for the appearance of a suppression in 
the DOS at the Fermi level. Previous wor k 10 i 43 was per- 
formed on a lattice appropriate for the description of 
LiAh/Ti2_ y 04, namely on a corner-sharing tetrahedral 
lattice. Also, that work was completed using a restricted 
HF theory: In terms of the above-presented formalism, 
the effective local fields hf were set to be zero, mean- 
ing that the spin degrees of freedom were restricted to 
develop along the quantization axis (z). 

All of our results here are for a 3d simple cubic lattice 
- unlike the corner-sharing tetrahedral lattice, this lat- 
tice is bipartite and therefore unfrustrated. In terms of 
the magnetic degrees of freedom, here we present results 
from our new HF calculations for two different situations. 
First, we further restrict the system to be paramagnetic, 
meaning that we do not allow for the formation of any 
local moments - this means that is forced to be equal 
to nil- Second, we relax the restriction of the moments 
being along the z-axis, and allow them to form in the 
a;z-plane, i.e., hf is nonzero but real. 

Some of our results for the DOS for a three-dimensional 
simple cubic lattice are shown in Fig. [TJ The lattice size 
is 18 3 , and the electronic filling factor is 1/4. The in- 
teraction strength is U/t = 11.5, where the DOS at the 
Fermi level has its largest suppression, and the disorder is 
modelled with a box distribution for a disorder strength 
of W/t = 6. The nonparamagnetic HF result shows a 
strong suppression of the DOS at the Fermi level; how- 



T 




(E-E p )/t 

FIG. 1: [Colour online] Density of states of a three- 
dimensional simple cubic lattice of size 18 3 for an electronic 
filling of 1/4. The interaction strength is U/t = 11.5, and 
the disorder strength is W/t — 6. The nonparamagnetic HF 
result (dashed curve) shows a strong suppression of the DOS 
at the Fermi level, whereas the paramagnetic HF result (solid 
curve) shows no such suppression. The broadening used cor- 
responds to d = 0.024t. 



ever, the paramagnetic HF result shows no suppression 
at all. 

Since we only found a suppression of the DOS with 
nonparamagnetic HF solutions, both with hf = (in 
RefiiS) and hf being real, the magnetic moments are 
essential to the suppression and to the metal-insulator 
transition (at least within the HF context). The moments 
found for the real hf HF ground state are strongly non- 
collinear - e.g., see our discussion in Ref. 30 . Therefore, 
the restriction of collinear moments employed in Rcfs. 10 
is not important to the suppression. Also, the type of 
lattice does not matter because the suppression is found 
both for the frustrated corner-sharing tetrahedral lattices 
(in Ref,—) and for the unfrustrated simple cubic lattices 
studied in this paper. Also, the type of disorder does not 
matter, since the suppression is found either with a quan- 
tum site percolation model (in Ref. 10 ) or with a uniform 
box distribution type of disorder. 

To better quantify the presence of the local moments 
we have calculated the average magnitude of the moment 
per electron. We define an Edwards-Anderson-like order 
parameter 

1 N " 1 
m=—^2\(Si)\, Si = -$ ia T aa >c ia < (13) 
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where Si is the spin operator on site i, f are the Pauli 
matrices, and N e (N s ) is the number of electrons (lat- 
tice sites). The quantity rn is similar to the Edwards- 
Anderson order parameter— in spin glass theory, which is 
used to distinguish between glass and nonglass phases, 45 
but here we use m as a characterization of local magnetic 
moment formation. 

We have considered many different parameter sets and 
lattice sites (see discussions below) , and our results for to 
are shown in Fig.O We see that arc no local moments in 
the noninteracting (U ft — 0) or weakly interacting sys- 
tems (U/t < 3) for W/t = 2 and 6. For the 18 3 lattice 
and U/t — 3 and W/t = 6, m is only about 0.0045, and 
for the 18 3 lattice and U/t = 3 and W/t = 2 we find 
to = 0.00002; therefore, this quantity is expected to be 
zero in the thermodynamic limit. For larger U/t we do 
find moments, and for a given interaction strength to in- 
creases as one increases the strength of the disorder. In 
all of the subsequent data that we show having a suppres- 
sion of the DOS at the Fermi level, and non-Drude-like 
conductivity, the average magnitude of the local mag- 
netic moments to is always nonzero. 
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FIG. 2: [Colour online] The average magnitude of the mag- 
netic moment per electron (m) vs. interaction strength (U/t) 
for a 12 3 lattice with disorder strength W/t = 6 (squares), and 
for an 18 3 lattice with W/t = 2 (diamonds) and 6 (circles). 



Now, we consider the variation of the DOS in more 
detail. We have calculated the DOS for a 1/4-filled three- 
dimensional simple cubic lattice of size 12 3 for various 
strengths of interaction U/t = 0, 1, 2, . . . , 8, 9, and 11.5. 
The on-site energies obey a uniform box distribution, and 
the disorder strength is W/t = 6. Results are averaged 
over 4 realizations of disorder - a limited study using 
more realizations of disorder found no qualitative changes 
when more realizations were employed. The maximum 
differences for the self-consistency are between 2 x 10~ D 
(for U/t = 1 and 2) and 1.2 x 10~ 3 (for U/t = 11.5). 
(Typically, the average differences were at least one order 
magnitude smaller than the maximum differences.) In 
Fig. [3l we show the DOS of this system obtained with d = 
0.048i. (The DOS curves are smooth and well separated 



from one another around the Fermi level, whereas for 
smaller broadening these curves are bumpy and not well 
separated, and we therefore choose d = 0.048i as the 
broadening width for the DOS of this system.) 

500 
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FIG. 3: [Colour online] Density of states of the Anderson- 
Hubbard model on a 1/ 4-filled three-dimensional simple cubic 
lattice of size 12 3 for a disorder strength of W/t = 6. Results 
without interactions (U/t = 0) and with interaction strengths 
U/t = 1, 2, 6, 9, and 11.5 are averaged over 4 realizations 
of spatial disorder. The vertical line represents the Fermi 
level. The broadening width is d — 0.048t. 



As shown in Fig. [31 for U/t = to 3, each DOS does not 
show any suppression. In fact, the DOS at the Fermi level 
is enhanced as the interaction is turned on and then in- 
creased. The suppression appears for the first time when 
U/t = 4, and the amount of the suppression increases as 
U/t is increased. We find that the maximum suppres- 
sion, for 1/4 filling, occurs around U/t = 11.5, although 
there is little difference between the DOS for this U/t 
and other Hubbard interactions close to this strength. 



C. Behaviour of the Optical Conductivity: 

Extracting the low- frequency behaviour of the conduc- 
tivity is nontrivial, and we briefly outline the numerical 
approach taken. We have examined the ac conductivity 
obtained with different values of the Gaussian broaden- 
ing width d in Eq. ([7]), and have selected a value for d 
based on the following: For a three-dimensional simple 
cubic lattice of size 12 3 with disorder strength W/t = 6 
and interaction strength U/t = 6, we calculated the ac 
conductivity with broadening widths d = 0.012£, 0.03i, 
and 0.06t (results not shown). Because the broadening 
width d has a finite value, the contributions of frequencies 
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between and d makes the imaginary part of the current- 
current correlation function ImA(q = 0, uS) = u)a\{uj) 
nonzero around lo — 0. As a result, when calculating 
<Ji(u>) by dividing ImA(q = 0, u>) by to, we obtain di- 
verging a\(ui) around u> = 0. The divergence is simply 
an artifact of numerical procedure used (a broadening 
width d that is finite), and is not associated with the 
physics of the system being investigated. Therefore, we 
have to cut off the conductivity curves at low frequencies 
(~ d) where they "turn up". Recalling that the DOS 
data of Fig. [3] used d = 0.048i, we choose d — 0.03£ as 
the broadening width for the conductivity because the 
corresponding conductivity curves are smooth and still 
retain much of the low-frequency behaviour - simply, we 
are not forced to discard as much low-frequency data as 
we would with d = 0.048£. (As reviewed in the discussion, 
this low-frequency behaviour is an important quantity to 
know.) 

We now discuss our results, juxtaposing DOS and ac 
conductivity data for each parameter set. First, we dis- 
cuss results for studies done on lattices of size 12 3 with a 
disorder strength of W/t — 6; an average over 4 realiza- 
tions of disorder is used for each data set. We increase 
the interaction strength U/t from zero up to roughly 12. 
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FIG. 4: [Colour online] Density of states (a) and ac conductiv- 
ity (b) for a 1/4- filled three-dimensional simple cubic lattice 
of size 12 3 with a disorder strength W/t = 6 with no inter- 
actions (U /t — 0) and with Hubbard interaction of strengths 
U/t = 1, 2, and 3. The DOS do not show any suppression at 
the Fermi level, and in fact increase as the interaction strength 
U/t increases. The low- frequency conductivity is enhanced as 
U/t is increased, and the ac conductivity is Drude-like. We 
have used a broadening corresponding to d = 0.0484 for the 
DOS and d = 0.03t for the conductivity. 



The results for the noninteracting case (U/t = 0) and 
for interactions U/t = 1, 2, and 3 are shown in Fig. 2] 
The solid lines in (a) and (b) represent the DOS and 
the ac conductivity for the noninteracting electrons, re- 
spectively. The ac conductivity has a shape that is typ- 
ical of a metal, but the low-frequency peak is broad. 
When we turn on the interaction to U/t = 1 (dotted 



lines), we see an enhancement in the DOS at the Fermi 
level. Concomitantly, the low-frequency ac conductiv- 
ity also increases from the noninteracting value. These 
enhancements could be a result of the screening of the 
disorder by the Hubbard interactions. As the interaction 
strength U/t increases to 2 (dashed lines) and 3 (dash- 
dotted lines), we find that the DOS at the Fermi level 
and the low-frequency ac conductivity are both further 
enhanced. Therefore, for a disorder of W/t — 6 at 1/4 
filling, a weak interaction of U/t < 3 enhances the low- 
frequency ac conductivity due to an increase in the DOS 
at the Fermi level. Recall, from Fig. [21 that for this range 
of interaction strengths no local moments are formed. 
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FIG. 5: [Colour online] Density of states (a) and ac conductiv- 
ity (b) for a 1/4- filled three-dimensional simple cubic lattice 
of size 12 3 with a disorder strength W/t — 6 and Hubbard 
interaction strengths U/t — 3, 4, 5, and 6. Each of the DOS 
for U/t = 4, 5, and 6 shows a suppression at the Fermi level, 
and the amount of the suppression increases as U /t increases. 
The low- frequency conductivity decreases as U/t increases, 
and the ac conductivity is no longer Drude-like. 



The results for interactions U/t = 3, 4, 5, and 6 are 
shown in Fig. [H A suppression of the DOS at the Fermi 
level first appears for U/t — 4, and its value is smaller 
than that of the noninteracting electrons and that of the 
U/t = 3 system. The low-frequency conductivity for 
U/t = 4 is also smaller than that of the U/t — 3 system. 
As the interaction strength U/t increases to the values 
of 5 and 6, the amount of the suppression of the DOS 
at the Fermi level increases, and the low-frequency con- 
ductivity decreases. Each of the ac conductivity curves 
still extrapolates to a nonzero value when the frequency 
u/t — > 0. We note that starting from U/t = 5, the low- 
frequency conductivity is smaller than that of the nonin- 
teracting electrons. More importantly, the ac conductiv- 
ity for U/t = 4, 5, and 6 is no longer Drude-like. In fact, 
our conductivity curves are similar to the localization- 
enhanced Drude theory of Mott and Kaveh,— and such 
conductivities have been observed experimentally previ- 
ously in the metallic phase near the metal-insulator tran- 
sition in experimental work on conducting polymers4£ 
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We show the results for interaction strengths U/t = 6, 
9, and 11.5 in Fig. © Compared to U/t = 6, the DOS 
at the Fermi level gets further suppressed. However, the 
amount of the suppression for U/t = 9 and 11.5 are not 
that different. It is very clear that, on the scale of Fig. [6l 
the dc conductivity for U/t = 9 and 11.5 are both zero. 
The ac conductivity for U/t = 9 and 11.5 both increase 
as uj/t increases up to the value of u/t ~ 3. 

We fit the ac conductivity for U/t = 6, 9, and 11.5 with 
power-law relations. For U/t = 6, the system appears 
to be metallic, and we fit its ac conductivity using the 
following equation: 

(7i H = we 2 t 2 (a + l3u a ), (14) 

where o~o, (3, and a are the parameters to be determined. 
The range of frequency over which we chose to fit the 
data is 0.04 < u/t < 1.2 (recall that we used d = 0.03< 
in producing these conductivity curves), and we obtain 

ai(u/t)/(TreH 2 ) = (0.002 ± 0.001) 
+ (0.145 ± 0.001)(u;/i) - 501±0 - 007 . (15) 

Here, the uncertainties are one standard deviation. We 
learn from this power-law fit that the dc conductivity for 
the system with W/t = U/t = 6 is <j = 0.002 ± 0.001, 
which is quite close to but still above zero. Therefore, we 
expect that the system is a metal. We also see that the 
exponent of the frequency is a — 0.501 ± 0.007, which 
is the exponent (0.5) that appears in the ac conductiv- 
ity of noninteracting electrons with strong disorder, in 
the metallic regime close to the transition^ where <ti(u) 
near u = goes as ai(u) — <jq ~ \fu. 

We fit the ac conductivity for U/t = 9 and 11.5 over 
the ranges of u/t e [0.074, 0.267] and u/t e [0.074,0.330], 
respectively. Setting cto = we obtain 

cn(u/t)/(TreH 2 ) = (0.23 ± 0.01) (a;/*) 1 ' 89 * - 02 (16) 

for U/t = 9 and 

ai{u/t)/{izeh 2 ) = (0.178 ± 0.006)(cjA) 2 - 05±0 - 02 (17) 

for U/t — 11.5. These exponents are very close to 2, 
and, again, this is the same as the exponent for nonin- 
teracting electrons for an even stronger disorder, namely 
for a disordered system in the insulating phased where 
<j\(u) ~ u 2 in the low- frequency regime. 

Therefore, the system with W/t = 6 is metallic for 
U/t < 6 and insulating for U -^9. By examining 
other U ft around 9 (data not shown) we can identify 
that the metal-to-insulator transition for this disorder 
and electronic filling indeed takes place at some critical 
(U/t) CT « 9. 

It is always desirable to complete studies with larger 
lattices, where finite-size effects are hopefully less puni- 
tive. In our study the benefits of having such data are 
(i) better energy resolution of the DOS near the Fermi 
level, and below we make clear the usefulness of such 
data; and (ii) better resolution of the ac conductivity 
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FIG. 6: [Colour online] Density of states (a) and ac conductiv- 
ity (b) for a 1/4- filled three-dimensional simple cubic lattice 
of size 12 3 with a disorder strength W/t = 6 and Hubbard 
interactions of strength U/t = 6, 9, and 11.5. Each DOS 
shows a suppression at the Fermi level, and the low-frequency 
conductivity is also suppressed - both suppressions increase 
with increasing interactions. The ac conductivity for U/t — 6 
extrapolates to a nonzero dc conductivity as w/( -> 0, how- 
ever, the dc conductivity for U/t = 9 and 11.5 are both very 
close to zero. 



at low frequencies. Using our computer resources the 
largest lattice that we have managed to treat, and are 
able to explore various parameters sets, is an 18 3 lattice 
with 1/4-filling. For this lattice with our HF formulation 
employing real hi, it takes about one month to achieve 
self-consistency in the HF calculations for a single re- 
alization of disorder, and due to this enormous amount 
of time we solve the problem only for one realization of 
disorder. 

For a disorder strength of W/t = 6 and interaction 
strengths of U/t = 0, 3, 6, and 9, we show the resulting 
DOS in Fig. [7] (Due to the higher number of energy 
levels per unit frequency we have used 1/2 the broadening 
width (d = 0.024i) as we did for the 12 3 studies.) We 
note that the behaviour of the DOS in this larger lattice 
is very similar to the 12 3 lattice that we have studied 
earlier. 

The corresponding ac conductivity for the same sys- 
tem is shown in Fig. [5] Note that this data corresponds 
to a very small Gaussian broadening of only d = O.Olt, 
and therefore we can obtain data down to very low fre- 
quencies. We clearly see that the data is similar (qualita- 
tively and quantitatively, indicating that our data does 
not suffer from strong finite-size effects) to that shown 
earlier for a 12 3 lattice. Further, the extrapolation of the 
conductivity in the zero frequency limit makes clear the 
nonconducting nature of the U/t = 9 system. 
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FIG. 7: [Colour online] Density of states for the Anderson- 
Hubbard model on a 1/ 4- filled three-dimensional simple cubic 
lattice of size 18 3 . For an intermediate disorder of W/t — 6, 
interaction strengths of U/t = 0, 3, 6, and 9 are considered. 
The Gaussian broadening width is d — 0.024£. 




FIG. 8: [Colour online] AC conductivity of the Anderson- 
Hubbard model on a 1/ 4- filled three-dimensional simple cubic 
lattice of size 18 3 for a disorder potential of strength W/t — 6, 
for interaction strengths of U/t = 0, 3, 6, and 9. The Gaussian 
broadening width is d — 0.0 It. 



D. System with a weak disorder: 

So far, we have studied systems with a disorder 
strength of W/t = 6. As mentioned earlier, this corre- 
sponds to a noninteracting disordered electronic system 
with a very low dc conductivity - that is, this corresponds 
to a "bad metal" . We found a suppression of the DOS at 
the Fermi level for interaction strengths U/t>4. When 
this happens, we find that the ac conductivity is no longer 
Drude-like, and for strong enough interactions, U/t ~ 9, 
the system no longer possesses a metallic conductivity. A 
natural question then arises: For smaller strengths of dis- 
order does one still obtain a metal-insulator transition? 

In this subsection we present our results for a weaker 



disorder corresponding to W/t — 2. Compared to 
the noninteracting W/t = 6 system, the noninteracting 
W/ 1 = 2 system has a mean free path and dc conductiv- 
ity almost an order of magnitude larger. 

For such disorder we are forced to use a large lattice 
because the energies are strongly degenerate in small lat- 
tices. In Fig. [5] we show the DOS for an 18 3 simple cubic 
lattice with interaction strengths U/t = 3, 6, 9, 11.5, 12, 
and 15. 
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FIG. 9: [Colour online] Density of states for a 1/4-filled three- 
dimensional simple cubic lattice of size 18 3 with a box distri- 
bution type of disorder of strength W/t = 2 (4 realizations of 
disorder) and Hubbard interaction of strengths U/t — 3, 6, 9, 
11.5, 12, and 15. 



For this weak disorder strength, the interaction 
strength of U/t — 3 does not produce a suppression in 
the DOS at the Fermi level. However, for U/t — 6, 9, 
11.5 and 12 a suppression in the DOS appears at the 
Fermi level, and the amount of suppression is deepest at 
U/t = 11.5. For U/t — 15 we find that the suppression 
has been eliminated. 

In Fig. [TU1 we show the ac conductivity for the same 
system with interaction strengths up to U/t = 12. The 
conductivity for U/t = 3 is clearly metallic and Drudc 
like, and extrapolates to roughly a^c ~ lAne 2 t 2 in the 
low- frequency limit. This value is roughly half that of 
noninteracting electrons with the same strength of disor- 
der. 

With increasing U/t the conductivity at low frequen- 
cies is lowered, with U/t = 9 and 11.5 having very sim- 
ilar behaviour. However, for U/t = 12 the conductiv- 
ity increases near u> = — this behaviour persists for 
even larger Hubbard energies. Therefore, this disorder 
strength does not lead to a nonmetallic conductivity for 
any Hubbard energy that we have studied. 

To aid in the understanding the above results, namely 
that for a larger disorder one does obtain a metal- 
insulator transition, whereas for smaller disorder one 
does not, we have examined the local charge densities 
for these parameter sets. In Fig. 111! we show histograms 
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FIG. 10: [Colour online] AC conductivity for a 1/4-fiUed 
three-dimensional simple cubic lattice of size 18 3 with a box 
distribution type of disorder of strength W/t — 2 (4 real- 
izations of disorder) and Hubbard interaction of strengths 
U/t = 3, 6, 9, 11.5, and 12. 



of the local charge density for the 18 3 lattice with disor- 
der strength W/t = 6 for U/t = and 9 (a), and W/t = 2 
for U/t = and 11.5 (b). 
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FIG. 11: [Colour online] Histograms of the local charge den- 
sity for the 18 3 lattice, (a) U/t = and 9 for a disorder 
W/t = 6. (b) U/t = and 11.5 for a weak disorder W/t = 2. 

For disorder strength W/t — 6, when the system is 
noninteracting the charge density is spread out over a 
large range, with one large peak close to n ~ 0.1. That 
is, in the absence of interactions and for this seemingly 



large disorder (see previous discussion) the ground state 
corresponds to a very small occupation of sites with large 
on-site energies, with a increasing occupation as the on- 
site energies decrease. However, when the interaction 
strength is strong (U/t = 9), the original peak at the 
small charge density does not change much, but a large 
new peak around n w 0.9 forms. That is, for this disorder 
the interactions lead to very different charge densities, 
and correspond to those of an insulating system. The 
contrasting situation for weak disorder is clear from the 
figure - the charge distribution is quite similar between 
the noninteracting (U/t = 0) and strongly interacting 
(U ft = 11.5) systems, and in both cases one has metallic 
conduction. 

We have also plotted the histograms for the effective 
on-site energies, that is ef^ from Eq. (TJ| averaged over 
both spins, since this quantity has been suggested to be 
important in the understanding of the metallicity of such 
ground statesi 49 i 50 However, the histograms for (ef^)a-, 
for the same parameter sets as those used in Fig. [TTJ 
show essentially no difference for these interacting HF 
ground states, one of which we found to be metallic and 
one insulating. 



IV. DISCUSSION 

Our numerical analysis of the real-space self-consistent 
HF treatment of the Anderson-Hubbard model away 
from 1/2 filling on 3d simple cubic lattices gives rise to 
three main conclusions, each of which we discuss in rela- 
tion to previously published theoretical and experimental 
work. 

• At least at 1/4 filling, and thereby well away from 
1/2 filling, one requires sufficiently strong disorder 
to obtain a metal-to-insulator transition. 

As preliminary work to our study, we first charac- 
terized the consequences of different strengths of disor- 
der for noninteracting electrons. For a disorder strength 
much less than the noninteracting bandwidth, given by 
B = I2t, we find a large dc conductivity, consistent with 
theoretical predictions of a large mean free path. How- 
ever, for a disorder strength that is 1/2 of this band- 
width, we find a greatly reduced dc conductivity, and 
applying Eq. (j 1 2|) , estimate a mean free path of the same 
order as the lattice constant. Therefore, for the latter dis- 
order strength the effects of the localization of electronic 
eigenstates is considerable. Note that we have treated 
the disorder of the Anderson-Hubbard model exactly, and 
have solved systems having large lattices with as many as 
18 3 sites, thereby allowing for such localization physics 
to be present in our HF ground states. 

At 1/4 filling and a disorder strength of W/t = 6, we 
find a critical interaction strength of U c /t » 9 leads to 
a metal-insulator transition, whereas a smaller disorder 
strength of W/t = 2 does not lead to any such transition 
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for any value of U jt that we have studied. Therefore, one 
concludes that a disorder that is sufficiently strong and 
thereby leads to strong localization effects is required to 
obtain such a transition. Of course, as always throughout 
our conclusions, this is what one may conclude via a HF 
treatment of this model. 

• In order to obtain a suppression of the DOS at the 
Fermi level, one must allow for the development of 
local magnetic moments. 

The importance of local moment formation in the dis- 
ordered metallic state was made clear in the seminal work 
of Milovanovic, Sachdev, and Bhatt^ The importance of 
such moments in obtaining the metal-insulator transition 
is demonstrated in our work. 

The results from Fig. Q] make clear the necessity of 
allowing for local moments (in a HF treatment) to form 
if one is to obtain a suppression of the DOS. Further, the 
results from Fig. [2 along with the DOS curves shown 
in Figs. HI [5] and El demonstrate that to obtain such a 
suppression of the DOS local magnetic moments must 
be present. The lattice and model of disorder do not 
influence whether or not such moments appear. 

We note that similar suppressions of the DOS within 
HF treatments of this and related models have been 
found by Logan and Tusch^ one of us and coworkers^ 
and most recently by Shinaoka and Imada^S In all of 
these studies a restricted HF formulation that required 
magnetic moments to form along some single chosen 
quantization axis was employed, again emphasizing that 
importance of allowing this degree of freedom. However, 
whether or not "twisted spins" that point in all directions 
possible are allowed^ , does not seem to be required. 

Due to the considerable interest in the DMFT treat- 
ment of correlated electrons, it seems appropriate to 
note that theories going beyond the single-site KKR-CPA 
treatments of disorder within DMFT, such as Refi^i, also 
find non-Fermi-liquid behaviour. More recently, a vari- 
ant of DMFT that, like this paper, also tries to treat 
the disorder exactly, and used the HF approximation in 
the evaluation of the off-diagonal self energy, has found 52 
zero bias anomalies similar to those shown in this paper. 

• When such moments develop, in addition to the 
suppression of the DOS, one obtains a non-Drude- 
like ac conductivity. With sufficiently strong dis- 
order and interactions, the dc conductivity is sup- 
pressed to zero and one obtains a metal-insulator 
transition. 

At least within our HF treatment, just above and just 
below the interaction strengths leading to the metal- 
insulator transition, one finds a — <7d c ~ V^J and a ~ to 2 
power law behaviour, the same as one finds for nonin- 
teracting electrons. Whether this is a by-product of HF 
being an effective one-electron theory remains as an out- 
standing question. 

We note that novel ac conductivities, including results 
that are beyond that which we conducted, namely that 



included the temperature dependencies, were also seen 
in the recent work of Kobayashi et a/. 35 These authors 
used the same HF decomposition as that used here. How- 
ever, unlike our work these authors focussed on 1/2-fillcd 
systems in two dimensions, and, in particular, on the 
properties of the novel 2d metallic phase that was pre- 
dicted based on earlier HF workj 32 i 33 i 34 These authors 
do not propose specific frequency dependencies, unlike 
the results we give in Eqs. I|14I16I17[) . so the potentially 
important effect of differing dimensionalities remains un- 
known. 

As mentioned in the introduction, the suppression of 
the DOS at the Fermi energy and the associated metal- 
insulator transition found in our earlier work on lithium 
titanate, 10 was the original motivation for this study. 
Since that paper new exact diagonalization and Monte 
Carlo results 2 ^ have made clear (albeit in two dimen- 
sions) that such physics is not an artifact of using a HF 
approximation. Subsequent HF work ) 20 ' 53 using averages 
over very large numbers of complexions of disorder, has 
found qualitatively similar results, but has also managed 
to produce many more energy eigenvalues close to the 
Fermi energy, and thereby better probe the functional 
form of the suppression of the DOS found in all of these 
papers. 

The origin of the suppression of the DOS is not deter- 
mined by the present study, although we have placed con- 
straints on what physics (local moments and sufficiently 
strong disorder) must be included in a model that will 
produce such behaviour. In terms of existing theories, 
we note that for weak disorder and weak interactions one 
may complete a derivation 5 ^ of the Altshuler and Aronov 
result^ 2 - for the Anderson-Hubbard model, and one finds 
that one gets an increase in the DOS as U/t increases 
from zero. That is, one finds 

""■' ■-WF '" r (18) 

where A is some positive constant and Nq(Ef) is the 
DOS of the noninteracting system at the Fermi energy. 
This is consistent with the W/t — 6 and U/t = 0,1,2 
and 3 DOS curves that we found for the 12 3 lattice, as 
well as the W/t = 6 and U/t — and 3 curves for the 18 3 
lattice (even though the weak disorder kf£ e ^> 1 assump- 
tion of the Altshuler and Aronov theory is not obeyed for 
our W/t = 6 system). For larger U/t our HF numerics 
show that the sign of the shift of the DOS changes, in 
that one has a suppression rather than enhancement, as 
is clear from Figs. 0] and O Given that we find that one 
must have local moments to have any suppression, even 
for weak disorder (see Figs. [2] and [9]), it's clear that one 
must extend the weak disorder paramagnetic Fermi Alt- 
shuler and Aronov theory, which remains an outstanding 
problem. 

However, a completely different idea was made in a 
recent HF and exact diagonalization stud y 20 ! 53 and these 
authors proposed that while the Altshuler and Aronov 
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form of 

N(E) ~\E-E F \ 5 (19) 

may be adequate away from the Fermi energy, very near 
the Fermi energy one instead finds a so-called soft Hub- 
bard gap. They showed that such physics in three di- 
mensions gives rise to 

N(E) ~ exp{-(- 7 log \E-Ep |) 3 } (20) 

For one-dimensional systems with strong localization ef- 
fects, these authors provide independent justification for 
their theory. However, to access the energies very close 
to Ep in both one and three dimensions, averages over 
very large numbers of complexions of disorder were car- 
ried out3i££ 

Our numerics do not have such good resolution, but we 
have examined DOS data for an 18 3 lattice for W/t = 6 
and U/t = 11.5, now averaged over three different com- 
plexions of disorder. First, we note that with a broad- 
ening roughly 1/5 of that shown in our Fig. [7] we find 
that the DOS indeed does indeed go to zero at Ep; this 
is consistent with Refsi 20 ' 53 . Second, fitting the DOS 
data with these smaller broadening factors, for ui/t > 0.2 
we find that Eq. (TT9|) fits the data well, and we obtain 
5 = 0.51 ±0.01. Therefore, we find the same exponent as 
in the theory of Altshuler and Aronov^ However, while 
our data for uj/t < 0.1 is sparse and therefore not com- 
pletely reliable, the functional form proposed in Eq. ([ST)]) 



indeed fits our DOS curves reasonably well. Our lim- 
ited numerics therefore support the soft Hubbard gap 
proposa l 20 ! 53 at energies very close to Ep. 

Finally, we mention the relationship of our work to 
some experimental results on transition metal oxides that 
were mentioned in the introduction. Sarma et al^i have 
investigated the electronic structure of LaNii-^Mn^Oa. 
The conductance is proportional to VV, where V is the 
bias voltage. As mentioned above, such a form fits our 
data very well. Kim et al^ have studied the transport 
and optical properties of SrTii-^Ru^Os. Their low fre- 
quency optical conductivity data is very similar to that 
which we obtained within HF. Our data also agrees with 
the optical conductivity found in the disordered metal- 
lic phase of conducting polymers SL The analysis of this 
data, based on the localization modified Drude modelp^ 
leads to the same 0.5 exponent as we obtained in our un- 
biased power law fitting, the result of which is given in 
Eq. (dU). 
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